Two-point correlation function in systems with van der Waals type interaction 



Daniel Dantchev 

Institut fur Theoretische Physik, Technische Hochschule Aachen, 52056 Aachen, Germany 

and 

Institute of Mechanics, Bulgarian Academy of Sciences, Acad. G. Bonchev St. Building 4, 1113 Sofia, Bulgaria 

(February 1, 2008) 

The behavior of the bulk two-point correlation function G(r;T\d) in d-dimensional system with 
van der Waals type interactions is investigated and its consequences on the finite-size scaling proper- 
ties of the susceptibility in such finite systems with periodic boundary conditions is discussed within 
mean-spherical model which is an example of Ornstein and Zernike type theory. The interaction 
is supposed to decay at large distances r as r~^ d+c '\ with 2 < d < 4, 2 < <r < 4 and d + a < 6. 
It is shown that G(r;T\d) decays as r - ^ -2 ' for 1 -C r <C £, exponentially for ^ < r <C r*, where 
r* = (cr — 2)£ln£, and again in a power law as r~^ d+a ^ for r 3> r*. The analytical form of the 
leading-order scaling function of G(r; T\d) in any of these regimes is derived. 
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I. INTRODUCTION 



It is well known that the critical properties of a given statistical-mechanical system depend only on a small number 
of parameters like the dimensionality d of the system, the symmetry of the order parameter characterizing the 
corresponding phase transition, and on general properties of the interaction coupling the order parameter at different 
locations. For example, in an isotropic 0(n) system one expects that all critical exponents and scaling function of a 
given physical quantity are independent on, say, lattice structure, or on short-range details of the interaction. Let us, 
for dcfmitcncss of notation, speak about Ising-like systems (i.e. n = 1) with Hamiltonian 

« = ~E J(*-*')SrSr>. (1.1) 

In the context of the critical phenomena the usual criterion for a given interaction to be considered as short-ranged is 
finite second moment of J(r), i.e. in terms of the Fourier transform J(k) of this interaction for small k — |k| one has 

J(k) = J + J 2 k 2 + A J(k), (1.2) 

where AJ(k) is asymptotically smaller than k 2 . Then, if 2 < d < 4, for the bulk two-point correlation function 

G(r; T\d) =:< S r S > - < S > 2 (1.3) 

one normally writes 

G(r; T\d) = D r - {d - 2+ ^ X ± {r / £), (1.4) 

where are two universal scaling functions (for T > T c and T < T c , respectively), D — D(T) is a nonuniversal, 
slowly varying function of temperature that well can be approximated by a constant in the vicinity of the critical 
point T — T c and £ is the bulk correlation length, i.e. £(T) ~ &t~ v , t -> 0+ with t = (T - T c )/T c . For T > T c one 
has 

X+a;(<i-3)/2+77exp(-a0, x -> oo 
X+(x) c { (1.5) 
constant, x — > 0. 

The above is, in fact, the classical result of Fisher Q for the two-point correlation function in the critical region of a 
simple fluid. Note, that when r 3> £ the correlations decay exponentially fast with the distance r. It is well known that 
in simple nonpolar fluids the interactions are characterized by potentials that decay as inverse powers of the distance 
at large r. In d = 3 the most prominent case is the induced dipole-induced dipole (or van der Waals) interaction for 
which (neglecting the retardation effects) the potential decays as $(r) := — J(r) = — A/r e , where A > is a positive 
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cons tant. One easily can check that this interaction has a Fourier transform which is indeed of the type given in Eq. 



(1.2). But for such systems (ci-dimensional Ising model), in which the interactions decrease in inverse power with 
respect to the distance between the interacting objects, i.e. as r ~^ d+a \ a > 2, the following rigorous result due to 
Iagolnitzer and Souillard || is available. 

Theorem j^]: The two-point correlation function G(r;T, H\d) of a ferromagnetic system in the presence of an 
external magnetic field H does not decay faster than its potential J. For any T < oo and any real H there exists a 
strictly positive constant C(T,H), such, that 

G(r; T, H\d) > C(T, H)J{r). (1.6) 

For H = the theorem is valid for T > T c and for any of the two "pure" phases (the "plus" and the "minus" ones) 
for T < T c . 



The immediate consequence of this theorem is that if T ^ T c (1.4) - (1.5) could no t be true for r large enough 



independently on how close T is t o th e critical point. The only way to reconcile (L4) - (L5) with the above theorem 



is to realise that if T ^ T c (1.4) - (1.5) could be valid only up to some r = r*(T). Then for 1 C r C ( the correlations 
will decay as r~( d-2+, '\ for £ <C r <C r* they will fall off exponentially, but, for r > r* they should again decay in a 
power law as a function of the distance, namely as r -( d + a ) . I n other words one should observe a crossover from power 
law to exponential and then, again, to power law behavior of the correlations. Saying this one immediately stacks 
with at least two important questions that appear naturally: 1) What is the value of r* , i.e. where this crossover 
happens and 2) What are the properties of the function describing that crossover. One of the aims of the current 
article is to answer those questions in the framework of an exactly solvable model. 

It is easy to check that one has the above situation only with interactions of the type J(r) ~ A/r d+(T , where a > 2. 
To avoid misunderstanding in the remainder of the text let us make the following definitions, i) An interaction will 
be called of short range if for any finite m its m-th moment is finite, i.e. '}2 ir r rn J(Y) < oo. ii) An interaction is long 
ranged if there exists a finite m such that the corresponding m-th moment diverges. If m = 2 this is a leading-order 
long range interaction, and if m > 2 this is a subleading (van der Waals type) long-range interaction. 

We recall that even if the interaction is short-ranged in the above sense (this is the situation we have with nearest 
neig hbo ur, next-nearest neighbour, etc. interactions, i.e. with interactions that are essentially of a finite range), then 



(1.4) - (1.5) are again valid only for 1 << r << r* r . The exact results for d = 2 Ising model (see, e.g. McCoy and 
Wu ||) and the mean-field results of Fisher and Burford Q suggest that, if T ^ T c , r* r ~ £ 2 . (The exact calculations 
due to Chen and Dohm j5| for the spherical model give a bit more "generous" estimation for r* r , namely r* r ~ £ 3 , 
see also below.) For r >> r* r the interactions decay, of course, again exponentially, but they contain a nouniversal 
prefactor ||], i.e. their leading-order behavior is then nonuniversal. If the interaction is of a leading long-range type 



then (1.4) is also valid but in the limit r/£ 3> 1 one has to require that X + (x) ~ X + x n ~ 2 ~ a = X + x~ 2cT , where we 
have taken into account that 77 = 2 — a if er<2 [|[j8| . The corrections to the large distance correlations in this case 
are in a power-law of r, which means that their leading order behavior is universal for any r >> 1. This asymptotic 
is confirmed by the exact results for the spherical model due to Joyce |||l(| and it is in tune with the above theorem 
for the Ising model. Note that for r large enough the correlations always fall off in a power law with the distance with 
the only exception of interactions of a fully finite range when they do decay exponentially. 



That (1.4) - ( p.. 5] ) should be modified for the case of subleading long-range interactions has been noticed by several 
authors. 

First Widom proposed [|ll| that for r — > 00 

G(r; T\d) ~ 0J(r) + a^-^' 2 exp(-r/a 2 ), (1.7) 

where (3 — l/(fegT), ks being the Boltzmann's constant, and ai and ai are "depending only on the thermodynamics 
state constants" . It is clear that in nowadays formulations the above means to take for the correlation function a sum of 



/3J(r) and the right-hand side of Eq. (1.4). Later the problem has been attacked by Enderby, Gaskell and March p2[. 
They consider a three-dimensional fluid, i.e. the case d = a = 3. Supposing the Ornstein-Zcrnikc integral equation to 
be valid and taking the direct correlation function to be c(r) = /3J(r), they obtain, after assuming that the structure 
factor S has a Fourier transform of the type S(k) = x/P + c 2 k 2 + c 3 k 3 + ■■-, that G(r; T|3) ~ J(r)x 2 /0 = A X 2 /(Pr 6 ), 
when r — > 00 and in temperature regions "well away" from the critical point J13|. Starting from this result, Kayser 



and Raveche |14| suggest that G(r;T\d) can be decomposed in two additive contributions G sr and G lr , where G sr is 
given by (1.4), plus higher-order terms that account for the usual corrections to scaling, and G lr = <d(r — r*)J(r)x 2 /P- 
Here Q(x) is the Heaviside step function and r* is to be determined by the requirement that at this point G sr = G . 
Taking \ ^ £ 2 ~ 1) in the expression for G lr , which in fact means supposing G(r;T\d) ~ J(r)% 2 //3, r — > 00, to be valid 
for general d and also for temperatures close to T c , the above authors derived 

r* = (a + In f (1.8) 
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In plj] Floter and Dietrich make similar statements for r* for the case d = a = 3. 

In the present article we will investigate the large r behavior of the correlations and will derive the explicit form of 
G lr within the mean spherical model. The interaction will be supposed to be of the type J = A/r d+a , with 2 < d < 4, 
2 < a < 4 and d + a < 6. 

If one knows G one can immediately determine the behavior of the bulk susceptibility \ by using the fluctuation- 
dissipation relationship x(T\d) = (3 ^ r G(r; T\d). Definitely, if the finite-size two-point correlation function is known 
for a given finite system with a characteristic size L, then one can determine in this way also the behavior of the 
finite-size susceptibility x{T\ L\d). In a recent article Q Chen and Dohm have addressed the question: could one say 
what should be the scaling structure of the finite-size susceptibility under periodic boundary conditions if one knows 
only the bulk two-point correlation function? They suggest a hypothesis, that this is possible by interpreting in a 
proper way the functional dependence of G on r as a dependence of x(T; L\d) on L. In the present article we check 
the relationship that they suggest between G and x(T; L\d) on the example of our exactly solvable model. For that 
aim we will use the results for x(T; L\d) derived in |17| for the same model. 

The structure of the article is as follows. In Section || we briefly describe the model and present our starting 
analytical expressions. Section [II contains our results for the large r behavior of the two-point correlation function 
G. In Appendix [A] we present some details of the calculations needed to determine the asymptotics of G as a function 
of r and £. Section |y| comments on the relationship between the derived results for G and the behavior of the finite- 
size susceptibility of systems with subleading long-range interactions. The article closes with a Discussion (Section 
0) where we speculate about the possible extensions of our results for other models. 



II. THE MODEL 



We consider a d-dimensional mean spherical model fig] , fl9| (for a comprehensive review on the results available 
for this model see pcjl ). The degr ees of freedom consist of a set of N localised spins with Gaussian weight, and the 



Hamiltonian is given by Eq. (IT). The interaction J(r) is supposed to be of van der Waals type, i.e. its Fourier 



transform is supposed to be of the form 

J(k) ~ J(0) (1 - v 2 k 2 + v a k° - v 4 k 4 + 0(k 6 )) , (2.1) 

where k = |k|, 4 > a > 2 and J(0), V2, v a and v 4 are nonuniversal positive constants. Note that the signs of the 
coefficients in the small k expansion of the Fourier transform of the interaction are chosen so as they normally appear 
for subleading long-range inte ractions that decay in power law with the distance between the interacting objects 
- molecules or spins. In ( [2.l| ) J(0), V2, v a and V4 are c-dependent — for simplicity of notation this dependence 
is omitted here. The term v a q a in ( |2.l| ) is associated with a contribution to the real-space interaction going as 
r -d-a F ur thermore, we suppose that J(k) — J(0) < if k ^ 0, which reflects the fact that there are no competing 
interactions in the system and that the only ground state is the ferromagnetic one. Of course, it would be interesting to 
consider such systems — say with a combination between antiferromagnetic short range and ferromagnetic subleading 
long-range interactions, but this is out of the scope of the current article. 
The partition function of the model is given by the multiple integral 



dsi--- I ds N exp[-/3H], (2.2) 

> 

supplemented by the mean spherical condition 



JV 

L 



^(4) = N, (2.3) 



which can be enforced with the use of a "Lagrange multiplier" term going as A s f t ne effective Hamiltonian, 
and thence into the partition function. The spherical model equation of state then takes the form 



, ^ = N. (2.4) 

A - J(0)(1 - v 2 k 2 + v a k° - v 4 k 4 ) 

The phase transition in this model occurs when the combination A — J takes on a value asymptotically close to zero. 



The difference between the equation of state in (2.4) and the standard mean spherical model condition in short range 
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systems lies in the addition of the term going as k a in the denominator on the left hand side of (2.4). In general, this 
term is taken to be negligible, but we will soon see that it leads to interesting effects. 

For the model defined in the above way it can be shown, following jl(|, that the bulk correlation function G(r; K\d, a) 
is given by, if 2 < d < 4, 



G(r;K\d, a) 



1 1 



r dr 



K (2ir) d J nd f- 2 + k 2 - bk a + ck A ' 



(2.5) 



where K = f3v 2 J(0), and b — v a /v2 > and c = V4/V2 > are nonuniversal constants. Let us note that the values 
of b and c are such, that there are no real roots of th e equation 1 — bk a ~ 2 + ck 2 = 0. The last follows from the 
propositions we made for J(k). Note also that in (2.5) we have taken the cut-off in the fc-space to be infinity (for a 
lattice system it will mean that one considers the limit of a zero lattice spacing). This is possible because of the rapid 
oscillations of the exponential function in the integrand, but in this way we neglect all finite cut-off effects that will 
give nonuniversal contributions towards the critical behavior of the two-point correlation function (see [|| for details). 
In (2.5) £ = £2 is the second moment correlation length defined via (see, e.g. EJ) 



G(0;K\d,a) 



^ ^G(k;K\d,a) 



(2.6) 



fe=0 



where G(k; K\d, a) is the Fourier transform of G(r; K \d, a). Because of this identification one can, in fact, skip for 
our purposes the analysis of the spherical field equation (2.4) - one directly has A = J(0)(1 + w 2 £~ 2 ). 



Since a > 2 and since we are i nter ested in the behavior of G(r; K\d, a) for |r| 1 (note t hat then the leading order 
contributions of the integral in ( |2.5P will be coming from small k values), one can rewrite ( |2.5| ) in the form 

(2.7) 



G(r; K\d, a) = G sr (r; K\d) + G lr (r; K\d, a), 



i.e. as a sum of "short-range" and "long-range" parts. The "short-range" correlation function is the part that is only 
due to the short-range component of the interaction and, as usual, will be taken to be of the form 



G sl \r;K\d) 



1 1 

K (2tt) c 



ik.r 



■R< 



k 2 ' 



(2.8) 



The other contributions that are due to the subleading components of the interaction do form the corresponding 
"long-range" part. As it has been already stated in Section | such a structure has been supposed to hold by Kayser 
and Raveche |l4| in their qualitative analysis of the correlation functions in fluids. 

Performing the integrations in (|2.8|) and taking into account that for 2 < a < 4 and r>l 



G(r:K\d,a) 



G sr (r;K\d) 
G sr {r-K\d) 



(&fc CT - ck 4 )e lk - r dr 



1 1 

K J2n) d J nd (£- 2 + k 2 - bk° + ck 4 )(£- 2 + k 2 ) 



k a 



r dr 



we obtain 



where 



G(r;K\d,a) 



K (2tt)<* J nd (£-2 + k 2 ) 2 

X sr {r/£) + br- {a -^X h (r/£_) + 



-(d-2) 



K (2TT) d / 2 



(2.9) 



(2.10) 



(2.11) 



sm[(d 



\x d l 2 +°-* [xl d/2 (x) -(d + a- 2)I d/2 _ 1 (x)] ) , 



(2.12) 



and • • • stays for contributions which are corrections with respect to the terms retained. Here I a (x) is the modified 
Bessel function, and p F q (a; b; z) is the regularized generalized hypergeometric function 
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P -F 9 (a;b; 



p F q (a;b;z) 



T(b 1 )T(b 2 )---T(b q y 
where p F q (&]b; z) is the generalized hypergeometric function 

(ai)fc(a 2 )k ■ • • (a p ) fe z k 



p F q (a;b;z) = ^ 



fe=0 



(6i)fc(62)fc---(6,)fc fc! 



(2.13) 



(2.14) 



The symbol (a)fc = a(a + 1) • • • (a + k — 1) = T(a + fc)/T(a) in the above equation is the Pochhammer's symbol. The 
function p F q is finite for all finite values of its arguments. In the above expressions only the leading order long-range 
contributions (i.e. the contributions "proportional to 6"), have been retained and we have supposed that 2 < d < 4, 
2 < cr < 4, and d + a < 6. 

We recall that for the Ornstein-Zernike type theorie s (inc luding the mean-spherical model, see, e.g. ]ic|] ) 77 = 0. In 
Section |y| we will di scuss briefly the generalization of ( 2.1C ) for models with 77 ^ 0. 

The expressions (pT0|)-(pT2|) are the analytical basis for our further analysis. Let us note that the correlations 
within the spherical model have been a subject of detailed investigations (see [[L0) and [^0| for a comprehensive review) 
for both short-range and leading long-range interactions. Surprisingly enough, they have never been investigated for 
subleading long-range interactions. 



III. LARGE DISTANCE BEHAVIOR OF THE BULK TWO-POINT CORRELATION FUNCTION 

The asymptotics of the scaling function X SI for 2 < d < 4 are well known (see, e.g. |2j| and references cited therein) 

' v /fa;( d - 3 )/ 2 exp(-a;)(l + 0(x- 1 )), x -> 00 
X sl '(x) ~ I (3.1) 
nd/2-i) zjWfzf + Otx 2 ) x^O 

2( 4 - d >/ 2 + 2 d / 2 sin(7rd/2)r(d/2) + ^ X X U ' 

Let us, nevertheless, make some comments here. First, let us note that the above asymptotic is obtained if one makes a 
quadratic approximation of the spectrum and lets the cut-off A of the theory go to infinity. (If one keeps a sharp finite 
cut-off A with such an approximation of the spectrum one will obtain a nonexponential oscillatory power-law behavior 
||, fl22|]). Second, on a lattice, for nearest neighbours interactions between the spins embedded in a d-dimensional 
cube it has been shown that the above expression is valid || only for 1 << r << £ 3 . If r > £ 3 the correlations 
do depend (up to the leading order) on the mutual positions of the spins involved, i.e. the lattice anisotropy comes 
into the play and can no longer be neglected [^J. So, one can think that the above expression is valid in the region 
1 << x << £ 2 . Third, for x — > the second term in this short-range expansion of the correlation function involves a 
t dependence of the type t a if one takes into account that £ ~ £q t~ v for t — > + and that (d — 2)v = a, i.e. 

*-fr/a - + + ° {i " ()2) ' ( * (3 ' 2) 

The temperature dependent term in this expansion is usually not explicitly specified in the literature on the spherical 
model. Finally, let us note that according to the above asymptotics and under the approximations made £2 = £e = 
i.e. the second-moment correlation length coincides in such a theory with the exponential-decay correlation length. 
(In H it has been shown that for a model on a hypercubic lattice £2 = a /[2 sinh(a/2£ e )], where a is the lattice spacing 
and £ e has been chosen to be along one of the principal axis of the lattice.) 

The asymptotics of X h for2<<i<4, 2<er<4, d + a < 6 are (see Appendix [X] for a derivation) 

f _ 2CT+rf /2-2 g (d +g -2)r((d+ g )/2-l) 3,-4 + 0(a .- 6)) 

X h (x) ~ { (3.3) 

0<r+d/2-4 r((d+ ff )/2-2) d+a-i 7r(rf/2+ CT /2+l) n(^\ t -» 

I, Z r(2-cr/2) x 2<*/ 3 r(d/2)sin(7r(d+<7)/2) ly V x ^ x ~ ' u - 

The above asymptotics lead to the following behavior of the short and long-range parts of the bulk two-point 
correlation function 



G*<(r;K\d) ~ --g^^-M/a^^-D/a exp (-r/0(l + O ((r/^)- 1 )) , r » (3.4) 
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and 



G h (r;K\d,a) 



b 2 



a-2 



K TT d / 2 



er(d 



r((d + a)/2-l) 4 (d+a) 

' r(i- ff /2) 4 ' >>4 ' 



(3.5) 
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FIG. 1. We present as illustration plots of the scaling functions X BT (x) and X lr (x) of the short-range and long-range 
correlation functions, respectively, as a function of the scaling variable x = r/f. X sr (x) is plotted for d — 2.5, whereas (x) is 
for d = 3 and a — 2.5. For any 2 < d < 4, X sr is always positive and decays monotonically as a function of x. For large values 
of x, X BT decays exponentially fast. Note that in contrast with X BT (x), X lr (x) is not a monotonic function of x. In addition 
X lr can be both positive and negative. X lT (x) decays in a power law, as x~ 4 , for large values of its argument. 

One can determine the crossover region where the correlations from short range become long range type. To that 
aim one has to solve the equation 



G SY {r-K\d) ~ G h (r;K\d,a). 
Having in mind Eqs. ( 2.10| ), ( |3.4| ) and ( |3.5[ ) one obtains that the crossover takes place at 



where 



a lnlnC 



(3.6) 



(3.7) 



The leading-order term of this result coincides with that one given in JT3] (if one takes into account that 77 = for 
the model under consideration). 

For d = a = 3, i.e. for the true van der Waals interaction, the corresponding scaling functions are 



(3.8) 



and 



G 



X h (x) = f~{l-~x [exp(-x)Ei(a;) - exp(a;)Ei(-x)] + ^x 2 [exp(-x)Ei(x) + exp(x)Ei(-a;)]| . (3.9) 
The asymptotics of the short-range scaling function are obvious, while these for the long-range one are 




(3.10) 



FIG. 2. The scaling function X lr (x) of the Ion g-range correlation function for d = a = 3. One observes that, in contrast with 
the short-range correlation function, it is not a monotonic function of the scaling variable x = r/£. It changes sign at x ~ 1.088 
and x ~ 6.146 and reaches a minimum at x m i n ~ 2.113 which is X r (a; m in) — —0.192. In other words the long range part of 
the interaction increases the correlations (in comparison with an effective short-range system having the same value of K; we 
recall that K is a a-dependent quantity) for r up to 1.088 £ and for r > 6.145 f, but decreases them for 1.088 £ < r < 6.146 £. 
The maximum of X lr is reached at x — and it is X lr (0) = y/2/ir ~ 0.798. The last implies, as it is to be expected, that for 
a fixed r the maximal increment of the correlations due to the long-range part of the interaction is reached at T = T c . X lr (x) 
decays in a power law, as x~ , for large values of its argument. 



The asymptotics of the correlation function at T = T c and for any hxed £ can be derived to much greater details 
in the limit r — > oo for this especially important case. They are (see Appendix ^) 



G(r;K\3,3) 



K it 2 r 6 



1 + 120(r/£)- 2 + 10080(1 - \p{r/S,)- A + 0((r/Z)- 6 ) 



(3.11) 



and 



G(r;K c |3,3) 



1 1 

K r Anr 



1 + — r- 1 - -b(b 2 - 2c)r- 3 + 0(r~ 5 ) 

TT 7T 



, r — ► oo. 



(3.12) 



One can easily check that up to the leading-order terms these asymptotics coincide with the corresponding ones that 
follow by using the behavior of the short- and long-r ange correlation functions given above. For example, when r — > oo 
but r < £ (i.e. x = r/£ -> ) from Eqs. (Oh and (O) one has 



G(r; if|3, 3) 



1 1 

if 47rr 



l + a; + ia; 2 + — r 1 - — r 1 x 2 bxx + 0(r~ 3 , x 3 , x 2 ^ 1 ) 

2 7T 7T 



oo,x = r/^0. (3.13) 



T he crossover from short-range to long-range type behavior happens at r — r* where r* , in full agreement with Eq. 
(3.7) for d = a = 3, is given by 



r* =Mln£ + 51nln£ + 



In In C 



(3.14) 
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Note that Eq. (3.7) was derived under the condition that d + a < 6 and its "analytical continuation" to d = a = 6 
is no t obv ious. It is nevertheless valid because for r >> £ the leading-ord er t erm in the behavior of G(r; K\3, 3) (see 
Eq. (3.11)) can be obtained from that one of G(r; K\d, a), given by Eq. (|3.5|), if one sets d = a = 3 there. 




FIG. 3. We give here a representative example of the total two-point correlation function G(r; K\d, a) as a function of r for 
d = a = 3 and K = b = 1, £ = 20. As it should be expected for a ferromagnetic system G > and decays monotonically as a 
function of r. 



IV. FINITE-SIZE SCALING SUSCEPTIBILITY 



One usually describes the critical behavior of finite systems in the framework of the finite-size scaling theory 
"!H III . The standard finite -size scaling is usually formulated in terms of only one reference length, namely the 



bulk correlation length £. The main statements of the theory are that 

i) The only relevant variable in terms of which the properties of the finite system depend in the neighbourhood of 
the bulk critical temperature T c is L/£. 

ii) The rounding of the phase transition in a given finite system sets in when L/£ = 0(1). 

The tacit assumption is that all other reference lengths (such as lattice spacings, inverse cut-off, etc.) will lead only 
to corrections in the above picture. In addition, by analogy with the bulk short-range systems it is supposed that if 
a > 2 the finite-size critical behavior will be that of the corresponding short-ranged finite-size systems (see, e.g. p8|), 
characterised by exponentially fast decay of the finite-size dependence of the thermodynamic quantities at least when 
the critical region of the system is leaved in the direction towards higher temperatures (the low-temperature behavior 
depends on additional features like existence, or not, of spin-wave excitations - Goldstone bosons). 

As it has been recently shown the above picture is, in fact, more complicated [p|]l7|,|29f| and not completely valid for 



systems with subleading long-range interactions |17|. For such systems within the mean-spherical model and under 
periodic boundary conditions it has been found that the finite-size susceptibility x(i; L) is of the form 

x(t;L) = L~</»Y(x l ,bL 2 -*), (4.1) 

or, equivalently, 

X (t; L) = U<l v [Y B *{ Xl ) + bL 2 -°Y h \ Xl )] , (4.2) 

where x\ = atL 1 ^ , and Y, Y SI and Y lr are universal functions. The quantities a, and b are nonuniversal constants. 
One quite common wa y of fi xing a is to choose it to be a = (Co") -1 ^- It is worthily to note the close similarity in 



the structures of Eqs. ( 2.1C ) and (4.2). In other words - if one knows the structure of the bulk two-point correlation 
function one easily can write the corresponding finite-size behavior of the susceptibility. A hypothesis about such a 
possibility has been stated for the first time in || . 

In the high-temperature, disordered phase, where tL x l v — > oo, we find that the long-range portion of the spin-spin 
interaction gives rise to contributions of the order of bL~( d+a ^ that swamp the exponentially small terms that are 
expected to characterise the signature of finite size in systems with periodic boundary conditions and short range 
interactions. In other words the subleading long-range part of the interaction gives rise to a dominant finite-size 
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dependence in this regime. This is entirely consistent with the inherent long-range correlations that accompany 
long-range interactions, but it violates the standard finite-size scaling formulation. More explicitly, one obtains 
Y SI (xi) ~ exp(— const, x"), while 

Y h ( Xl )~x^ 2 \ (4.3) 

when x\ — > oo. This asymptotic follows from the requirement the finite-size corrections to be of the order of L _ ( rf + cr ) in 
this regime, which is to be expected on general grounds and is supported by the existing both exact and perturbative 



results for models with leading long-range interaction included [|30|-|32|. Note that (L3) implies for the temperature 
dependence of this correction 

x{t\L)- X {t\oo)~t- dv -^L- {d+a \ tL 1/l/ ->oo. (4.4) 

Obviously, the existence of such power-law finite-size dependent dominant terms above T c is of significance in the 
analysis of Monte Carlo data for such systems. 

Let us now consider the case d + a — 6, 2 < d < 4, 2 < a < 4, which contains the genuine van der Waals interaction 



with d = a = 3. Then, instead of (4.2), one has [p/fl 



X (t; L) = L?/" {y»(xx) + bL 2 -° [Y*( Xl ) ln(i) + Y?( Xl )] } . (4.5) 

Comparing with the corresponding results for the correlation function one observes, since there is no explicit In r 
dependence there, that this subtle feature like the logarithmic-in-L corrections will not be captured in the above 
mentioned approach — from the structure of the bulk two-point correlation function with respect to r to obtain 
that one of the finite size susceptibility with respect to L (simply by considering L in the role of r). Nevertheless, 
indications that the situation here may be more complicated are found in the short distance expansion of the bulk 
correlation function which has logarithmic in r terms. 



V. CONCLUDING REMARKS AND DISCUSSION 



In the present article we derived the analytical behavior of the two-point correlation function in a system with van 
der Waals type interaction. The treatment has been made within the mean-spherical model, which is an example of 
Ornstcin and Zcrnike type theory. We have pointed out that the leading order behavior of G(r; T\d) as a function 
of the distance is exponential only within the region of separations r between the interacting objects given by the 
condition £ -c r <C r* = (a — 2)£ln£. Obviously, taking into account the dependence of £ on the temperature, this 
region widens essentially only very close to T = T c . When r is outside the region defined above the correlations decay 
in a power law as a function of r: as r~( d ~ 2 ) for r -C £, and as r~( d+<7 ^ for r 3> (a — 2)£ln£. It turns out that 
G(r;T\d) can be decomposed in a "short-range" and "long-ra nge" parts (see Eq. (p/7|)). The corresponding short- 
and long-range scaling functions X SI and X lr are given in Eqs. ( 2.11 ) and ( 2.12| ), respectively. For the case 2 < d < 4, 
2 < a < 4 and d + a < 6 the behavior of these function s is illustrated in Fig. whereas the small and large value 
asymptotics of the functions are given in Eqs. ( |3.l| ) and (3^). A special attention is paid to the most important case 
of d = a = 3 which mimics the real van der Waals interaction in fluids. The analytical expressions for the scaling 
functions are given in Eqs. (3.8) and (3.9). The behavior of X lr is plotted on Fig. | The asymptotics of G(r; K\3, 3) 
and G(r; K c \3,3) are given in Eqs. (3.11) and (3.12). The behavior of the total correlation function G(r; K\3, 3) is 
illustrated in Fig. |[ 

Let us note that since all of the above results are for the mean sph erical mode l th ey pertain to the case of r\ = 
models. Naturally, one stacks with the question: How expressions like (2.10) and (3.7) should be modified for models 
with 7/^0? A hint in this direction can be found in [ 
terminology, r* = (a — 2 + ry)£ln£. In order to reconci 
r\ ^ (see also Eq, 



14[ ] - for such models Kayser and Raveche suggest that, in our 
c this statement with Eq. (2.10) one has to suppose that when 



G(r;K\d,a) = D(T)r 



-(d-2+77) 



X 



db,sr 



(r/£) + br- {a - 2+ ^X ± ' h \r/C) + 



(5.1) 



Here A sr is supposed to have the usual properties (see Eq. (1.5)), whereas for X lr we suppose that X lr (x) — > X 







and X h (x) — > X+x 2 ' 2 *>', x 00, where X^l and are positive constants. The large value asymptotics of X h (x) 
ensures that the correlation function decays as r~( d+cr ) for r ^> r* , which is in full agreement with fll.6|), and that 
r* = (a — 2 + ?7)£ ln£, which coincides with the result of Kayser and Raveche |l4j. The property G(r; T\d) ~ J(r)x 2 / (3, 
v — > 00 lllfl is retained too. 
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We emphasize, nevertheless, that despite all of the above features, for the moment (5.1) is only a plausible hypothesis 
the verification of which is still lacking. 

At the very end we note that, according to a recent hypothesis j5j, the behavior of the bulk two-point correlation 
function G(r; K\d, a) can be related to that one of the finite-size susceptibility under periodic boundary conditions. 
In the present article (see Section IV) we checked this hypothesis within the mean spherical model using the results 
for the finite-size susceptibility derived in |i7|| . Definitely, one can extend the calculations presented here to models 
with 77 7^ by using renormalization group techniques. 
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APPENDIX A: MATHEMATICAL DETAILS 

In this Section we will provide the mathematical details needed to derive Eqs. (^3), ( 3.1C ), ( |3.1l| ) and ( 3.12 ). 
Let us start with the case d = a = 3. Then 



r)dk, 



where 



f(k) 



k 



£- 2 + k 2 - bk 3 + ck 4 ' 



(Al) 



(A2) 



It is easy to show that if for a given integer n > the derivatives f( p \Q) of the function / do exist for p = 0, • • • ,4n + 2 
and, in addition, f^(oo) = 0, p = 0, • • • , 4n + 2, then 



f(k)sm(kr)dk = r' 1 ^ [/ (4p) (0)r 

p=0 



-ip _ ^(4p+2) f n } r -4p-2 



(0)r- 



0(r 



-4p-5 



(A3) 



Applying this to (Al) in the li mit one immediately obtains (3.11). 

In order to derive Eq. (3.12) let us note that (£~ 2 = at K = K^) 



G(r;K c \3,3) = — 



111 
1 1 1 



1 



sin(fcr) 



dk 



l-bk + ck 2 k 

/>oo 

fb{k) sm(kr)dk — c j f c (k) sin(kr)dk 
o 



(A4) 



where 



f b (k) :-- 



l-bk + ck 2 ' 



(A5) 



and 



fc(k) := 



l-bk + ck 2 '' 



(A6) 



and use has been made of the fact that J °° sin (fc)/fc = 7r/2. Applying again (A3) for the evaluation of the integrals 
in Eq. (A4), we obtain t he res ult given in Eq. ( 3.12 ). 

The derivation of Eq. ( fTl~C|) is a bit more complicated First, let us note that 



G(r;if|3,3) = ^G(r; K c \3, 3) 



K 



G(r;X c |3,3) 



£- 2 + k 2 - bk 3 + ck 4 k 2 - bk 3 + ck 4 
1 



1 1 1 

K2n 2 r J Q (1 - bk + ck 2 )(£- 2 + k 2 - bk 3 + ck 4 ) k 



sin(fcr) 
k 

sin(fcr" 



dk 



dk 



^G(r;K c \3,3) 



k 3 



±_L!l 

K2it 2 r 
sin(Ax) 



sin(fcr) 



o (1 - bk + ck 2 )(£- 2 + k 2 ) 



dk- 



{£- 2 + k 2 ) 2 k 



di- 



ll 



§ G ^|3, 3) 



sin(fer) 



o r 2 + fc 2 



dk + b 



o t- 2 + k 2 



ain.(kr)dk- 



k 2 



sin(fcr) 



dk 

(C- 2 +P) 2 k 



K, 



o ^ 2 + fc 2 



d£- 2 Jo S,- 2 + k 2 



sm(kr)dk 



(A7) 



Above we have already dealt with the large distance asymptotic of G(r; iT c |3, 3). In order to obtain ( 3. 10|) now it only 
remains to note that Ip3 



1 sin(fcr) 7r 
p + fc2 jfc dfc = 2 



1 — exp — 



and 



/■°° sin(fc) 1 

/ 2 / 2 ^fe = tt - [exp(-x)Ei(x) - exp(x)Ei(-x)] , Re(x) > 0. 



(A8) 



(A9) 



At the end, let us derive the results given in Eq. ( |3.3| ). The case 1 <C r <C £, i.e. x — > 0, is simple - usin g the s eries 
representations of the modified Bessel functions I a (x) (see, e.g. Q or p4j| ) one obtains from Eqs. ( 2.12 )-( 2.14 ) the 
asymptotics of X ll (x) for small values of the argument, given in Eq. (3.3). Much more interesting is the case when 
r » ( > 1, i.e. when We will present here a derivation of the corresponding asymptotic of X h '(x) without 

making use of the large value asymptotic results for the function 1F2. One can get an impression of the beauty of the 
proposed way of acting only after taking a look at the results available for the function 1F2 (see, e.g., p5[ ). We start 
by noting that 



poo ,a+d/2 



d f c 



dt 



j-cr+d/2 
X + t 2 



Jd/2-l(t) 



= (l+x— ) / dz exp(-zx) / dt exp(-zt 2 )t ff+d/2 " 2 J d / 2 -i(t) 
ox J Jo 

In order to evaluate the last integral in ( AlOj ) one can use the formula |34j 

[3 v T{v/2 + n/2 + 1/2) 



dx x A 'exp(— ax )J„(/3x) = 
With its help one obtains 



; u+1: 

2 '4a 



r ( d+a — 'i) 



d + a d 1 
~~ 2 ; 2'~4z 



(A10) 



(All) 



(A12) 



r(^-i) 



2 d/2 T (|) 



-^(1 + x— ) / exp(-zx)z 1 - (d+CT)/2 iFi 

" ox 1 



1; 2 ; 4z 



(A13) 



Note now that when x ^> 1 the main contribution of the integral in the above expression will stem from small z 
values. Using the corresponding asymptotic [p3| 



L Fi(a;6; -y) 



T(b - a) 



y- a (l + a ( a -6+l) 2/ - 1 + 0(y- 2 )) 



(A14) 



of i-Fl(a; b, —y) for y ^ 1 and performing the integrations, we arrive at the asymptotic of A lr (x) reported in Eq. ( |3.3| ) 
for the case i>1. 
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